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Abstract 

We study Radiative Electroweak Symmetry Breaking in the Minimal Supersymmetric 
Standard Model (MSSM). We employ the 2- loop Renormalization Group equations for 
running masses and couplings taking into account sparticle threshold effects. The de- 
coupling of each particle below its threshold is realized by a step function in all one-loop 
Renormalization Group equations (RGE). This program requires the calculation of all 
wavefunction, vertex and mass renormalizations for all particles involved. Adapting our 
numerical routines to take care of the succesive decoupling of each particle below its 
threshold, we compute the mass spectrum of sparticles and Higgses consistent with the 
existing experimental constraints. The effect of the threshold corrections is in general 
of the same order of magnitude as the two-loop contributions with the exception of the 
heavy Higgses and those neutralino and chargino states that are nearly Higgsinos for 
large values of the parameter \i. 
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1. Introduction 



The low energy values of the three gauge coupling constants known to the present 
experimental accuracy rule out the simplest versions of the Grand Unified Theories. 
In contrast, supersymmetric unification, in the framework of Supersymmetric Grand 
Unified TheoriesQ @, is in excellent agreement @ with a unification energy scale Mqut 
within the proton decay lower bounds. Moreover, softly broken super symmetry, possibly 
resulting from an underlying Superstring framework, could lead to SU(2) L x U(l)y 
gauge symmetry breaking through radiative corrections for a certain range of values of 
the existing free parametersfl. In such a scenario the elegant ideas of Supersymmetry, 
Unification and Radiative Symmetry Breaking are realized within the same framework. 
The Minimal Supersymmetric extension of the Standard Model (MSSM)i incorporates 
all the above. Due to its minimal content and the radiatively induced symmetry breaking 
it is the most predictive of analogous theories. 

As in the numerous@0 existing analyses of radiative symmetry breaking in the 
MSSM, in the present article we employ the Renormalization Group. The Higgs boson 
running mass-squared matrix, although positive definite at large energy scales of the 
order of Mqut, yields a negative eigenvalue at low energies causing the spontaneous 
breakdown of the electroweak symmetry. The "running" of mass parameters from large to 
low scales is equivalent to computing leading logarithmic radiative corrections. Although 
this scenario depends on the values of few (3 or 4) free parameters, one could interpret 
it as leading to the prediction of Mz in terms of Mgut-, ot the Planck mass, and the 
top quark Yukawa coupling. Another way to interpret the predictions of this model is to 
consider Mz determined in terms of the supersymmetry breaking scale. The analysis of 
the results helps us find out to which extent the low energy data can constrain the type 
and scale of supersymmetry breaking. 
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The purpose of the present article is to include in the above stated scenario the 
so-called low-energy "threshold effects". Since we have employed the DR scheme in 
writting down the one-loop Renormalization Group equations, which is by definition 
mass-indepedent, we could "run" them from M GUT down to M z without taking notice 
of the numerous sparticle thresholds existing in the neighborhood of the supersymmetry 
breaking scale near and above M z . This approach of working in the "full" theory con- 
sisting of particles with masses varying over 1-2 orders of magnitude has to overcome the 
technical problems of the determination of the pole masses. Our approach, also shared 
by other analyses, is to introduce a succession of effective theories defined as the theories 
resulting after we functionally integrate out all heavy degrees of freedom at each particle 
threshold. Above and below each physical threshold we write down the Renormalization 
Group equations in the DR scheme only with the degrees of freedom that are light in 
each case. This is realized by the use of a theta function at each physical threshold. The 
integration of the Renormalization Group equations in the "step approximation" keeps 
the logarithms ln(^) and neglects constant terms. The physical masses are determined 
by the condition m(m p hys) = rn v hys which coincides with the pole condition if we keep 
leading logarithms and neglect constant terms. The great advantage of this approach is 
that the last step of determining the physical mass presents no extra technical problem 
and it is trivially incorporated in the integration of the Renormalization Group equa- 
tions. 

2. The softly broken Minimal Supersymmetric Standard Model 

The superpotential of the minimal supersymmetric extension of the standard model, 
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or just MSSM, is 

W = Y e L j E c H\e ij + Y d Q^D c a E\e l3 + Y u Q> a U c a H^ + I'll' 1 1 U,, (1) 

(e 12 = +1) in terms of the quark Q(3,2,l/6), D c {3, 1, 1/3), U c {3, 1, -2/3), lepton 
L(l,2,-l/2), £ c (l,l,l) and Higgs ^(1,2,-1/2), # 2 (1, 2, 1/2) chiral superfields. We 
have suppressed family indices. The second Higgs doublet H 2 is necessary in order 
to give mass to the up quarks since the conjugate of Hi cannot be used due to the 
analyticity of the superpotential. It is also required in order to cancel the new anomalies 
generated by the fermions in H\. Note that the superpotential (1) is not the most general 
SU(3)c x SU(2)l x U(1)y -invariant superpotential that can be written in terms of the 
given chiral superfields since terms like D C D C U C ,QLD C ,... etc not containing any ordinary 
standard model interaction could be present. The superpotential (1) could be arrived 
at by a straightforword supersymmetrization of the standard three Yukawa interaction 
terms. It possesses an anomalous R-parity broken by supersymmetry-breaking gaugino 
masses down to a discrete R-parity which ascribes -1 to matter and +1 to Higgses. There 
is also an unwanted continuous PQ-type symmetry leading to an observable electroweak 
axion which is broken by the last term in (1). This term introduces a scale fi which has to 
be of the order of the soft supersymmetry breaking scale in order to achieve electroweak 
breaking at the observed Mz value. Although it appears unnatural that the scale of 
breaking of a PQ symmetry should be related to the supersymmetry breaking, there exist 
schemes based on an enlarged framework (extra fields or non-minimal supergravitational 
couplings) that lead to dynamical explanation of the order ot magnitude of the scale ,J . 

The fact that supersymmetry is not observed at low energies requires the introduc- 
tion of extra supersymmetry breaking interactions. This is achieved by adding to the 
Lagrangian density, defined by the given SU(3)c x SU{2) L x U(l)y gauge symmetry 
and W, extra interaction terms that respect the gauge symmetry but break supersym- 
metry. This breaking however should be such that no quadratic divergences appear and 
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the technical "solution" to the hierarchy problem is not spoiled. Such terms are gener- 
ally termed "soft". The most general supersymmetry breaking interaction Lagrangian 
resulting from spontaneously broken Supergravity in the flat limit (Mp^oo) contains 
just four types of soft terms , i.e gaugino masses, <3>*<3>-scalar masses, $$$-scalar cubic 
superpotential interactions and $$-scalar quadratic superpotential interactions. For the 
MSSM this amounts to 

£sb = -^^AXA-m^lH^-m^^-mllQ^-mllDf - m 2 & \U c \ 2 

1 A 

- m\\L\ 2 - m%\Ef - (Y e A e U E c E\e l3 + Y d A d & a D c a E[e tJ 

+ Y u A u Q ja U c a W 2 e t3 + h.c) - (B^iJjif^ + h.c) (2) 

Again we have suppressed family indices. We denote with Ei,E 2 the ordinary Higgs 
boson doublets and with Q,D C ,U C ,L,E C the squark and slepton scalar fields .The gaug- 
inos \a are considered as four component Majorana spinors. Apart from the three 
gaugino masses, B and the two soft Higgs masses, there are also 5N G masses and 3N G 
dimensionfull cubic couplings in the simplest case that we retain only family diagonal 
couplings. Thus, totally 6 + 8N G new parameters. Note that these new parameters are 
dimensionfull and that without a simplifying principle they could in general represent 
different scales. 

A dramatic simplification of the structure of the supersymmetry breaking interactions 
is provided either by Grand Unification assumptions or by Superstrings. For example, 
577(5) unification implies at tree level mQ = rrifj — m^, mi — , M 1 = M 2 = M 3 and 
Ad = A e . SO (10) unification implies further equality of all sparticle masses, equality of 
Higgs masses and equality of the three types of cubic couplings. The simplest possible 
choice at tree level is to take all sparticle and Higgs masses equal to a common mass 
parameter m , all gaugino masses equal to some parameter m\/2 and all cubic couplings 
flavour blind and equal to A . This situation is common in the effective Supergravity 
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theories resulting from Superstrings but there exist more complicated alternatives. For 
example Superstrings with massless string modes of different modular weights lead to 
different sparticle masses at tree level!!. The equality of gaugino masses can also be 
circumvented in an effective supergravity theory with a suitable non-minimal gauge ki- 
netic term@. Note however that such non-minimal alternatives like flavour dependent 
sparticle masses are constrained by limits on FCNC processes. In what follows we shall 
consider this simplest case of four parameters m , m 1 / 2 , A Q and B Q . 



3. Radiative corrections and symmetry breaking 



The scalar potential of the model is a sum of three terms, being the contribution of 
F-terms, 

+ Y,\ Y dH{^D c a + Y n HiejiUZ\ 2 + El^ey^' | 2 + £|y u #^%f (3) 
where the sums are over the omitted family indices, the contribution of the D-terms 
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and finally the scalar part of —C-sb shown in (2). The tree level scalar potential leads to 



electroweak breaking as long as 

m\m\ - fi 2 B 2 < 

This is clear from 

V = m^H^ 2 + ml\H 2 \ 2 + (B^H[Hie ij + h.c) 

+ l -g ,2 {m 2 - \H 2 \ 2 ) 2 + l -g\\H x \* + |# 2 | 4 + A\HlH 2 \ 2 - 2|tf 1 | 2 |tf 2 | 2 ) + (5) 



written in terms of 



m? 2 = mL+ /i 2 - (6) 



We have replaced the appearing parameters with their running values mf(Q), mf>((5),... 
as defined by the Renormalization Group. The results based on the tree level study of 
the potential cannot be always trusted since they are sensitive to the choice of the renor- 



malization scale. Adding the one-loop radiative corrections obtained in the DR scheme, 

AVl = ^Str{.M 4 (M-M7Q 2 ) - 3/2)} (7) 

we end up with an Effective Potential that upon minimization supports a vacuum with 
spontaneously broken electroweak symmetry^"— . A reasonable approximation to (7) 



would be to allow only for the dominant top-stop loops. Note that although the Renor- 
malization Group improved tree level potential depends on the scale Q this is not the 
case for the full 1-loop Effective Potential which is Q-independent up to, irrelevant for 
minimization, Q-dependent but field-independent terms. 

Minimization of the 1-loop Effective Potential gives two conditions 

1 2 m\ - m 2 tan 2 (5 

2 Mz ~ tan* (3 - 1 (8) 



and 



sin2(3 = - _? BfJ _ 9 (9) 
mf + m,2 



The angle f3 is defined as (3 = tan 1 (i^/fi) in terms of the Higgs v.e.v.'s 
v\ y 2 =< H° 2 >■ The masses appearing in (8) and (9) are defined as 



m?, 2 = ml, + (10) 



All parameters are Q-dependent. At Q = Mz 

Ml= l -tf + f){vl + vl) + ... (11) 

Mz denotes the Z - boson pole mass Mz = 91.187CreV, and the ellipses are higher order 
corrections. Note also that in our convention fiB has the same sign with m\ + m\. 

We shall assume that at a very high energy scale Mqut the soft supersymmetry 
breaking is represented by four parameters m , m^, A Q and B of which we shall consider 
as input parameters only the first three and treat BiMz) as determined through equation 
(9). Actually we can treat j3(Mz) as input parameter and both B(Mz), fi(Mz) are 
determined by solving the minimization conditions (8) and (9), with the sign of \x left 



undetermined. The top-quark massO, or equivalently the top-quark Yukawa coupling, 
although localized in a small range of values should also be considered as an input 
parameter since the sparticle spectrum and the occurrance of symmetry breaking itself 
is sensitive to its value. Thus, the input parameters are m , m^, A a , (3(M Z ) and m t (M z ). 

Since radiative corrections are generally expected to be small with the exception of 
the contributions from the top-stop system, a reasonable approximation of (7) is 

AVx = E rntiHml/Q 2 ) - 3/2] - ^ 4[|n( 2 , g2) _ 3/2] (12) 



where we have kept only the t,t c and t contributions The field-dependent 
"masses" appearing in (12) are 

m t = Y t H° 2 , ml = ^{ml L + ml R ±[(ml L -m 2 RR ) 2 + Am A RL ^} (13) 
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where 

m\ L = m 2 i + ml + {^g' 2 -\g 2 )m\ 2 -\H° x \ 2 ) 

ml R = ml + ml- l -g>\\H°\ 2 -\Htf) 

m\ R = -Y t (H°A t + fMH?) (14) 

m RL = rn 2 LR* (15) 

Note that the use of the t-t contribution can be misleading in some cases due to large 
cancellations occuring with terms that are not included [Arnowitt and Nath in ref.6].A 
complete analysis requires that the contributions of all sectors to the effective potential 
are duly taken into account. 



4. The Renormalization Group and Threshold effects 

Consider the Renormalization Group equation for a soft mass parameter derived in 
the DR schemeEH 

dm b 

= ^ m (16) 
This equation should be integrated from a superlarge scale Q = M GUT , where we impose 

a boundary condition m(McuT) — m o, down to any desirable value of Q. As we come 

down from Mqut as long as we are at scales larger than the heaviest particle in the 

spectrum we include in b contributions from all the particles in the MSSM. When we 

cross the heaviest particle threshold we switch and compute b in a new theory, an effective 

field theoryEH with the heaviest particle integrated out. Coming further down in energy 

we encounter the next particle threshold at which point we switch again to a new effective 

field theory with the two heaviest particles integrated out. It is clear how we proceed 

from then on. 
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The change in the running mass parameter m at a particle threshold M in the above 
scheme comes out to be for m < M 

*H K %-^ n ,^) (17) 

where h + and 6_ are the Renormalization Group coefficients computed in the effective 
theories above and below the threshold respectively. Comparing (17) with the exact 
result obtained from the 2-point function associated with m we find that there is a finite 
non- logarithmic part that is missed by our approximation". The further M and m 



are apart the better the approximation becomes. Of course, the great advantage of the 
approximation lies in the fact that it is done entirely at the level of the Renormalization 
Group without the need to calculate the finite parts of n-point functions. 

The Renormalization Group equation (16) referring to a particular running mass 
m(Q) is integrated stepwise in the above stated manner down to the physical mass 
corresponding to m(Q). The physical mass is determined by the condition 

m(m phys ) = m phys (18) 



Note that in the DR scheme the inverse two-point function corresponding to the running 
mass m(Q) will be at 1-loop of the general form 



Y^\Q\ m(Q)) = Q 2 ( Cl + c 2 Zn(^^)) + m 2 {Q){c\ + c>(^))) (19) 



Imposing the condition (18) we see that the right hand side of (19) gives (ci + d^)rr? vhya . 
Thus, condition (18) coincides with the true (pole) condition for the physical mass only 
when the constant non- logarithmic contributions can be neglected. 

In what follows we shall present the 1-loop /3-functions of gauge and Yukawa couplings 
as well as those for the soft masses, cubic parameters A and quadratic parameters B 
and fJcM. Note that the threshold corrections introduced in our approximation by the 
theta-functions at 1-loop are expected to be comparable to the standard 2-loop RG 



corrections. In our numerical analysis that will follow we shall employ the 2-loop RG 
equations which have not been presented here due to their complicated form but can be 



found elsewhereOjJ. In our notation, for a physical mass M, 

9 M = 9(Q 2 - M 2 ) (20) 

Also t stands for t = InQ 2 and (3\ = ^ for each parameter A. Note also that we assume 
diagonal couplings in family space. 

The /3-functions for the three gauge couplings are 

§ s m = ^f T ' 3 ' • i = w (21) 

The coefficients b{ are 3 respectively. 

Keeping the Yukawa couplings Yt,b,r of the third generation fermions, the correspond- 
ing (3 functions are, 

= f 3 (Y T ) = ^{--T T2 g 2 2 --T Tl g l 2 + 2T TT Y T 2 + -Y b 2 } (22) 
dt v ' (4tt) 2 2 y 10 2 5 y ' 

dY b 



dt 



P(Y b ) = 

} b ;{~T b3 g 3 2 ~ lr b2 g 2 2 - ^-T bl9l 2 + \T bt Y t 2 + 3T bb Y b 2 + ]-Y T 2 } (23) 



(4tt) 2l 3 2"*"' 30 



f s W) = ^{-|W-|W-^W + sr.it> + iT.n'} (24) 

The threshold coefficients Tj, T T i,etc appearing in the expressions above are shown in 
Table I. We denote by G,W, B the S77(3),S77(2) and £7(1) gauge fermions respectively. 
The /^-functions for the cubic couplings are 

^ = ^{-3^ 2 2 M 2 ^ i -^ 1 2 M 1 (2 + ^ i )^ 

+ 3Y b 2 A b e 5 Q + AY T 2 A T + A T [Z Tigi 2 + Z T2 g 2 2 + Z TT Y T 2 ]} (25) 
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+ Y T 2 A T 6z~ L + Y t 2 AtO H2 u + 6Y b 2 A b 

+ A h [Z m g 3 2 + Z b2 g 2 2 + Z bl9l 2 + Z bt Y t 2 + Z bb Y b 2 }} (26) 
+ A[Z i3 (73 2 + ^2 2 + ^i 2 + ^ 2 + ^r5 2 ]} (27) 



In our notation 6 ab = 6 a 6 b and 9 abc = 9 a 9 b 9 c . The coefficients Z qi are shown in Table II. 

Next we procceed to the RG equations for the scalar masses. The RG equations for 

the sparticle masses refer to the third generation. For the other two generations the 

Yukawa couplings could be set to zero due to their smallness. 

din 2 - 1 & *\ 1 

= 7^{-[^3 2 (^-^) + ^ 2 2 (^-M + ^i%-^)]^| 

- ™gs 2 MSe 6 - 3g 2 2 M 2 2 9 w - ^ 9l 2 M 2 9^ + ± 9l 2 S 

+ Y 2 [m\B^ + ml9 D + m 2 2 9 H2 + A 2 9 H20 + ^{B Hxtj - 29^)}} 

+ Y b 2 [m%9 ki + m %6 5 + m\9 Hl + A 2 b 9 Hi3 + /i 2 ^ - 29^)]} (28) 
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^{-[^s 2 (% - e e ) + ^i 2 (% - 9b)] m 2 
f 9s 2 Ml9 d - f^ 2 Mle B - 2 - gi *S 

2Y t \m%d^ + + m 2 0„ 2 + A\Q H ^ + - 2^)]} 

^{-[^s 2 (^ - ^) + ^<7i 2 (*5 - 9b)\ m% 
™g 3 W 3 6 - ^M\B B + ^i 2 S 

2^ 2 K% + mj0Q + m 2 ^ + ^ fll<5 + /i 2 (^ 2(5 - 2^)]} 
^{-[^ 2 2 (% - 6 W ) + A^ 2 ^ - 6> 5 )] m| 
S^Mf^ - ^x 2 M 2 ^ - 

Y T 2 [m 2 L e 6l + + ml6 Hl + A 2 T 6 Hi ^ + // 2 (0^ - 2^)]} 

fsi 2 M 2 ^ + lg?S 



+ 2Y T *[m%e 6l + m\9 L + m\6 Hl + A 2 T 9 HiL + /z 2 (^ z - 26^)]} (32) 



dm{ 1 3 2 3 



{-[^ 2 2 (^ - + TH^ 2 (^ " ^b)]™? 



dt (4tt) 2L l 2^ v JJ1 10 



- 3^ 2 2 (M 2 2 + /i 2 )%^ - + - ^9i 2 S 



+ Y T 2 [m 2 1 + ml9 L + ml9^ + A 2 T 9^} 

+ 3y b 2 K + m|^ + m^ + A 2 ^]+3y t V%i/} (33) 
- 3£ 2 2 (M 2 + i?)Q h ^ - \g?{M\ + f)e 62§ + 



+ 3F t 2 [m 2 2 + + m 2 ^ + A 2 ^] 

+ SnV^ + nV^z) (34) 

The quantity S appearing in the equations above is defined as 

S= Tr{^e m m 2 } (35) 

In the absence of the threshold effects this quantity is multiplicatively renormalized. 
Therefore if it vanishes at the unification scale, due to appropriate boundary conditions, 
it vanishes everywhere and its effect can be ommited altogether from the RGE's. However 
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in our case this does not any longer hold owing to its explicit threshold dependence and 
S starts becoming nonvanishing as soon as we pass the heaviest of the thresholds. For 
the Higgs and Higgsino mixing parameters m\ = B\i and \i respectively we have, 

3 

+ ^ 2 + \v h * + \YA ml 



+ 3A t F t 2 ^ + 3A b n 2 ^ + A T F T 2 ^]} (36) 
3 2 



+ (^Hi + d H 2 ~ 8 ®HiH 2 + 0tfiB + ^ 2 s) 



^ 2 (^ + ^) + ^ 2 (^ + %) + ^ 2 (% + ^)}/^ (37) 



Finally the beta functions for the three gaugino masses are 

where 6j are the beta function coefficients of the gauge couplings given earlier and Si are 
threshold function coefficients given by, 

1 N 9 

Ss = -3^--E(%+^, + y (39) 
D i=i 
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1 1 

& = -6^-^E( 3 % + ^)-^(^+^ 2 ) (40) 
z i=i z 

= n[E(^ + ^ + ^ + ^ + %) + X^ + ^ )] (41) 

The dimensionful parameters, masses and cubic couplings, are meant to freeze out when 
the energy crosses below the mass scale associated with the heaviest particle participat- 
ing. This can be implemented by multiplying the corresponding quantity by the relevant 
theta function. Thus for instance A t freezes out below the thresholds of either t, or t c , 
or if 2, whichever is the heaviest, and the associated theta functions should multiply the 
r.h.s of eq. (27). For simplicity of our notation we do not indicate that explicitly in the 
RGE's displayed in eqs. (25)-(38). 

5. Formulation of the problem and numerical analysis 

The problem at hand consists in finding the physical masses of the presently un- 
observed particles, i.e. squarks, sleptons, Higgses, Higgsinos and gauginos, as well as 
their physical couplings to other observed particles. This will be achieved by integrating 
the Renormalization Group equations from a superheavy scale Mqut, taken to be in 
the neighbourhood of 10 16 GeV, down to a scale Q Q in the stepwise manner stated. If 
the equation at hand is the Renormalization Group equation for a particular running 
mass m(Q), then Q Q is the corresponding physical mass determined by the condition 
m (Qo) = Qo- If the equation at hand is the Renormalization Group equation for a gauge 
or Yukawa coupling the integration will be continued down to Q Q = Mz- Acceptable 
solutions should satisfy the constraints (8) and (9) at M z , i.e. describe a low energy 
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theory with broken electroweak symmetry at the right value of M z — 91.187GeV". 

The boundary condition at high energy will be chosen as simple as possible, postpon- 
ing for elsewhere the study of more complicated alternatives. Thus at the (unification) 
point M GUT , taken to be 10 16 GeV, we shall take 

mQ(M GUT ) = m^c{M GUT ) = m 0c (M GUT ) = m~ L (M GUT ) = m^ c (M GUT ) 

= m Hl (M GUT ) = m H2 (M GUT ) = m (42) 

and 

M^Mgut) = M 2 (M GUT ) = M 3 (M GUT ) = mi/2 (43) 
In addition we take equal cubic couplings at M GUT , i.e. 

A e (M GUT ) = A d (M GUT ) = A U (M GUT ) = A (44) 

All our boundary conditions are family blind. We shall also denote with B and \x Q the 
boundary values at M G \jt of the parameters B(Q) and n{Q)- This five parameters m , 
m i/2, A Q , B D and \i Q are not all free due to conditions (8) and (9) which could be viewed 
as determining B and as trading /i for f3 = tan^^/vi). Thus, m Q , m^, A a , (3{M Z ) as 
well as the sign of //(Mz) can be our free parameters. 

Our set of constraints includes the low energy experimental gauge coupling values 
which we have taken to be a 3 (M z )jfg = 0.117 ± 0.010, a em {M z Y x ^jg = 127.9 ± 0.1 
and (sin 2 6 w ) | w = 0.2316 - .88 10- 7 (M t 2 - 160 2 )GeV- 2 ^. The WS values for the 
couplings are related to their URQones through the relations gj^ = —Cg 2 /967r 2 ), 
where C = 0,2,3 respectively for the three gauge groups. The unification scale M G ut 
is determined from the intersection of olyjjr and gauge couplings and is found to 

be in the vicinity of 10 16 GeV. This value of M G ut is not easily reconcilable with the 

low energy value of a 3 quoted above and the universal boundary condition, as given in 

1 Note that at the 2-loop order the D R scheme needs to be modified so that no contribution to the 
scalar masses due to the "e-scalars" El shows up. 
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eqs. (42-44) if the effective SUSY breaking scale Ms is below 1 Tevl^l If we treat the 
low energy value of a.3 as an output we find it to be >0.125, i.e slightly larger than 
the most favourite experimental value quoted previously with a tendency to decrease as 
Ms gets larger than >1 Tev.Our interest in this paper is mainly focused on the mass 
spectrum and on the effect of the mass thresholds to it. The subtle issue of the gauge 
coupling unification in conjuction with the small value of shall be addressed to in 
a forthcoming publication. In the course of our numerical computations we allow for 
switching off the effect of the the thresholds from masses and cubic couplings involved. 
In this way we can compare the predictions for couplings and pole masses in the two 
cases i) With all thresholds present in both couplings and dimensionful parameters and 
ii) thresholds appearing only in gauge and Yukawa couplings. The latter case has already 
been considered by several groups. Although the difference is expected to be small only 
case (i) represents a consistent prediction of the spectrum, together with olgvt and 
Mqut, in the framework of the leading logarithmic approximation. 

For the b quark and r lepton our inputs are M£ ole = 5 GeV and M^ ole = 1.777 GeV 



which are related to their running MS masses by 

ro ^> = (1 + s (45) 

m r (M r ) = M T (46) 

These are evolved up to M z according to the SU C (3) x U em (l) RGE's and are then 



converted to DR values. From these the DR values for the Yukawa couplings Yb(Mz), 
Y T {Mz) are determined. The numerical procedure for the determination of Mqut re- 
specting the experimental inputs for a em ,sin 2 6w and M^ e needs several iterations to 
reach convergence. 



As for the top Yukawa coupling our input is the DR running mass m t {Mz) related 
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to the physical top mass Mf° 6 by 



mt(M t ) 



M t 



(47) 



The recent evidence 1111 for the top-quark 



mass has motivated values for M t in the neigh- 



borhood of 176 ± 8 GeV. 

In our numerical procedure we follow a two loop renormalization group analysis for 
all parameters involved, i.e couplings and dimensionful parameters, v.e.v's included. We 
start with the MS values for the gauge couplings at Mz giving a trial input value for 
the strong coupling constant «3 in the vicinity of .120, which are then converted to their 
DR values. These are run down to .\//,. M with the SU C (3) x U em (l) RGE' s to know 
the running bottom and tau masses as given before. From these by use of the RGE's 
for the running masses m&, m T we run upwards to know their MS values at Mz which 
are subsequently converted to their corresponding DR values. This provides us with the 
bottom and tau Yukawa couplings at the scale Mz- The top Yukawa coupling is known 
from the input running top quark mass m t (Mz) as said earlier. The evolution of all 
couplings from Mz running upwards to high energies determines the unification scale 
Mqut and the value of the unification coupling acur by 



Running down from Mqut to Mz the trial input value for o>3 has now changed. This 
procedure is iterated several times until convergence is reached. In each iteration the 
values of B, fi, which as stated previously are not inputs in this approach, are determined 
by minimizing the scalar potential. For their determination at the scale Mz we take into 
account the one loop corrected scalar potential. This procedure modifies the tree level 
values B(Mz), ^l(Mz)- It is well known that the value of /i affects the predictions for 
the physical masses especially those of the neutralinos and charginos. In approaches in 
which the effect of the thresholds is ignored in the RGE's the determination of B,fi is 



aiDji(M GUT ) = a 2 - m {M GUT ) = a GUT . 



(48) 
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greatly facilitated by the near decoupling of these parameters from the rest of the RGE's. 
However with the effects of the thresholds taken into account such a decoupling no longer 
holds since the thresholds themselves depend on -B,/i, or equivalently on fi, mf. 

In solving the system of the 33 RGE's involved, among these those for the v.e.v's, 
we have used special FORTRAN routines of the IMSL library available to us which 
use the Runge — Kutta — Verner sixth order method and are capable of handling stiff 
systems of nonlinear differential equations with high accuracy. 

Throughout our analysis we avoid considering values for tan (3 for which the couplings 
are driven to large values outside of the perturbative regime. 

The experimental lower bounds for the masses of new particles extracted from ac- 
celerator data are as follows^. The four neutralino mass eigenstates have to be heavier 
than 20,45,70 and 108 GeV while the two chargino states have lower mass limits of 43 
and 99 GeV. Charged sleptons have to be heavier than 45 GeV, while sneutrinos have to 
be heavier than 41.8 GeV. There is also a 150 GeV lower bound on the mass of squarks 
and gluinos. Charged Higgses should be heavier than 41 GeV while the CP-odd neutral 
Higgs should be heavier than 22 GeV. The lightest of the two neutral CP-even Higgs 
eigenstates should have a mass larger than 44 if the CP-odd Higgs is lighter than Mz or 
60 GeV if the opposite is true respectively. 

The interesting part of the output of the numerical integration of the RG equations 
consists of the mass spectrum of the new particles (gaugino-Higgsinos, squarks and 
sleptons) as well as the Higgs masses. The neutralino mass eigenstates can be read 



off from reference [JI2" 1. The two stop mass eigenstates t\ and t 2 correspond to the mass 
eigenvalues m| shown in(13).The Higgs mass eigenvalues at tree level are m\ = m\ + m\ 
for the neutral pseudoscalar Higgs A, and 

m 2 H , h = \[M z 2 + m 2 A ± y/(m\ + Mf ) 2 - 4Mfm> S 2 (2/?)] (49) 



2 We have also constrained the output to cases of neutral and colourless LSP. 
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for the two neutral scalar Higgses h and H. The charged Higgs has a tree level mass 
Mfj+ = m 2 A + Myy- As is well known, the 1-loop radiative corrections, mostly due to the 
large value of the top Yukawa coupling, are important for the Higgs masses. Following 



an approximation that has been tested in the literature 119, we have computed the Higgs 
mass eigenvalues based on the 1-loop effective potential (12), where the dominant third 
generation contribution has been kept. 

The radiative corrections to the pseudoscalar and charged Higgses are known to have 
an explicit dependence on the scale Q which cancels against the implicit Q dependence 
of the mixing mass m|(Q) and sin2/3(Q). Following other authorsB we choose as 
appropriate scale for the evaluation of their masses a scale Qq for which the effect of 
the radiative corrections is vanishingly small. Then Q turns out to be in the vicinity 
of the heavier of the stops and in this regime the stops have not been decoupled yet. It 
is therefore permissible to consider their loop effects to the effective potential as given 
in the references cited above. Keeping only the dominant third generation corrections 
to the Higgs masses the aforementioned cancellation among the Q dependent pieces is 
rather incomplete since the evolution of m\(Q) depends on the gaugino masses as well 
(see eq.(36)). Therefore there is a residual scale dependence, induced by the gauginos, 
whose effect on mf(Q) is not guaranteed to be small especially for large values of the 
soft mass m^. 

In our approach we have verified that the radiatively corrected Higgs masses for the 
heavy Higgses as these are calculated at the scale Mz and the corresponding masses at 
Qo are very close to each other even for large values of the soft mass mi/2- This is due 
to the appearance of the thresholds within the RGE of m\(Q) which properly takes care 
of the gaugino decoupling. This would not have been the case in schemes in which such 
a decoupling is not present. In those cases it is required that either the physical Higgs 
masses are evaluated as poles of the one loop propagators or the effect of the gaugino 
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fields is duly taken into account in the effective potential approach, for the effect of the 
radiative corrections to be numerically insensitive to the choice of the scale. This subtle 
issue is under investigation and the results of this analysis will appear in a forthcoming 
publication® . 



6. Conclusions 

We have displayed some of our results in tables III to VIII. In all tables we have 
taken fi(M z ) positive. The negative fi(M z ) values lead to qualitatively similar results. 
In the first four tables (III - VI) for a fixed value of m t (Mz) = 175 GeV, we present 
acceptable spectra that have been obtained for different values of the four parameter 
m , A Q , mi/2 and tanf3(Mz)- In the table III, for a characteristic set of values A Q = 400 
GeV, m = 300 GeV and m^j = 200 GeV we have varied tanf3 between 2 and 25. For 
values of tanfi larger than about 60 no electroweak breaking occurs in this case. Note 
the well knownEH approximate equality between the masses of one of the neutralinos 
and one of the charginos. The lightest Higgs turns out to be heavier than the Z - boson 
and increases with increasing the value of tan/3. Althought not displayed, for negative 
/i its mass drops below Mz for small values of the angle tan(3 ~ 2. In table IV, for 
fixed characteristic values of tanft = 10, A Q = 400 GeV and m = 300GeV,we vary 
mi/2 between 75 GeV and 700 GeV. The sensitivity of the whole spectrum on mi/2 is 
apparent. In table V, for fixed tan[3 = 10, m = 300 GeV and mi/ 2 = 200 GeV, we vary 
A Q between and 800 GeV. In this case, as in all other acceptable cases, the LSP is a 
neutralino with a mass roughly independent of A Q . We have also included in this table 
negative values of A in the range -800 GeV to 0. Note that in both cases the gluino mass 
is roughly stable. In table VI we keep tan(3 = 10, A Q = 400GeV and = 200GeV 
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fixed and vary m between 100 and 800 GeV. In all cases shown the LSP is a neutralino. 
In table VII we have presented spectra obtained, for fixed m , mi/2, A Q and tanf3 values, 
when we vary the top quark mass. In all cases displayed in the tables III - VIII we have 



the approximate equality m x o ~ m x cll£) due to the fact that fi(M z ) turns out to be 



2 A-2 

significantly heavier than Mz- 

As is apparent from the results displayed in tables IV to VI the value of the strong 
coupling constant has the tendency to decrease with increasing the supersymmetry break- 
ing scale. In the table IV for instance keeping fixed the values of m , A , a 3 (M z ) takes 
values between ~ .131 and ~ .127 if m^j varies from 75 to 700GeV. The situation is 
similar for the other two cases of tables V and VI where two of the soft SUSY breaking 
parameters are kept fixed and the third varies from low to large values. One notices 
however that the decrease in (Mz) in these cases is rather slower as compared to that 
of table IV. Thus (Mz) decreases faster in the direction of increasing the soft gaugino 
mass mi/2. 

Finally in table VIII for a characteristic choice of the parameter values we compare 
outputs for three dinstict cases. One loop predictions (case [a]), two loop predictions with 
thresholds only in couplings (case [b]), and the complete two loop case where thresholds 
appear in the RGE's of both couplings and dimensionful quantities (case [c]). Comparing 
cases [a] and [b] we see that as expected the value of the strong coupling is increased from 
its one loop value from ~ .118 to .133. The bulk of this increase, as is well known, is due 
to the two loop effects. At the same time the value of the unification scale is increased 
too from 2.1881 10 16 to 2.8876 10 ie GeV. Two loop effects delay the merging of the 
oi\ and a 2 gauge couplings resulting in larger values of a^. As can be seen from the 
comparison of the two outputs, two loops and threshold effects in couplings have a small 
effect on the mass spectrum. Differences are small of the order of 2% or smaller in 
all sectors except the neutralino and chargino states X121 X21 whose masses are exactly 
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Mi, M2 in the absence of electroweak symmetry breaking effects. In these states we have 
differences of the order of 10% or so originating mainly from the different evolutions of 
the soft gaugino masses Mi j2 and the 011,2 gauge couplings in the I - loop and 2 - loop 
cases. 

Switching on the the mass and cubic coupling thresholds produces a minor effect, as 
can be seen by comparison of cases [b] and [c], except in the case of the heavy Higgses 
and those neutralino and chargino states that are nearly Higgsinos, i.e. of mass ~ //, 
for large values of the parameter \i. These states are labelled as Xg 4 , xl I* 1 those cases 
relatively large differences are observed. For the neutralinos and charginos the large 
differences are attributed to the value of the parameter \x{Mz) which turns out to be 
larger, by about 10%, in case [c] as compared to that of case [b]. As far as the heavy 
Higgses are concerned from the discussion in the previous section it is evident that this 
discrepancy is due mainly to the evolution of m§, whose value affects substantially the 
masses of the pseudoscalar and charged Higgses, and in particular on its dependence 
on the gaugino masses. In order to determine the Higgs masses we need, among other 
things, the value of m\ at a scale Q a , which minimizes the loop corrections due to the 
quarks and squarks of the third family, knowing its value at M Z - However its evolution 
in the two approaches, with and without thresholds in masses and cubic couplings, is 
quite different. In the complete 2 - loop case, which takes into account all threshold 
effects, we have properly taken care of the decoupling of all particles as we move below 
their thresholds, gauginos included, unlike the case where the threshold effects are totally 
ignored. This results in rather large logarithmic corrections and hence to the relatively 
large differences (~ 6%) occuring in this cases. For a proper treatment of the radiative 
corrections in the second scheme, where threshold effects are not present in the RGE's, 
explicit one loop calculations of the gaugino contributions to the two point Green's func- 
tions of the Higgses should be carried out. It is only in this case where comparison of 
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the predictions for the heavy Higgses spectrum of the two approaches at hand can be 
made. At any rate the rather large differences seen in some particular cases point to the 
fact that a more refined analysis of the radiative effects to the Higgs sector is needed 
which also takes into account of the contributions of the gauginos and not just those of 
the heavy quarks. We have undertaken such a calculation and the results will appear in 
a future publication cQL 
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Table Captions 

Table I: Threshold coefficients appearing in the renormalization group equations 
of the gauge and Yukawa couplings. Above all thresholds these become equal to unity. 

Table II: Threshold coefficients appearing in the renormalization group equations 
of the trilinear scalar couplings. Above all thresholds these are vanishing. 

Table III: MSSM predictions for m t = 175 GeV, A Q = 400 GeV,m = 300, m 1/2 = 
200GeV and for values of tan/3 ranging from 2 to 25. Only the \i > case is displayed. 
The MS values of a em , sin 2 9w, a 3 and the unification scale are also shown. ctcuT is the 
DR value of the unification coupling. M t is the physical top quark mass, g denotes the 
gluino and xl, Xi are neutralino and chargino states. (ti i2 , 61,2) and (t 12 , v t ) are squarks 
and sleptons of the third generation while Ui,u\, di,d\ and ej,e£, Vi denote squarks and 
sleptons of the first two generations. h ,H , A, H ± are the CP even, CP odd and the 
charged Higgses respectively. 

Table IV: Same as in table III for m t = 175 GeV, A = 400 GeV, m = 300, tan/3 = 10 
and for values of m^j ranging from 75GeV^ to 700GeV (// > 0). 

Table V: Same as in table III for m t = 175 GeV, m = 300 GeV, m 1/2 = 200GeV 
and for values of A Q 0, ±300, ±800 GeV (/i > 0). 

Table VI: Same as in table III for m t = 175 GeV, A Q = 400 GeV, m 1/2 = 200 GeV, tan (3 = 
10 and values of m Q from lOOGeV to 800GeV (// > 0). 
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Table VII: Mass spectrum of the MSSM for input values m = 400 GeV, mi/ 2 = 
300GeV, A Q = 200 GeV, tan (3 = 20 and values for the running top quark mass equal to 
175, 180, 185, 190 and 195 GeV (/i > 0) respectively. Couplings and masses shown are as 
in table III. 

Table VIII: MSSM mass spectrum for the inputs shown in the first row (/i > 0). We 
compare 1 - loop (case [a]), 2 - loop with thresholds in couplings (case [b]) and com- 
plete 2 - loop predictions (case [c]) with thresholds in both couplings and dimensionful 
parameters. 



31 



TABLE I 



Ti = ^[20 + 6 ki + 6fj 2 + \{6 Hl +0 H2 ) + Y?i=i(\Q~ Li + + \^Qi + i°u t + \ & b z )\ 

T 3 = I - l G ~ i^ELll 2 ^ + o&i + 

T t2 = \[-l + 46» Hl - - 6>£^ + ^Oh^lw] 

T t i = ^[11 - 46>be + 86*^^ - 26^ + 46^ - - 46» Mi5i ] 

T tt = | [2 + O^e + 36 Hl + 26» Z ^J 

Tb3 = I [6 - 0&D - 6^] 

^62 = + 40 Hl - 20^^ - 6q W + ^h^vA 

T bl = ±[-21 - 40^ - 180^ + 240^ + 360^ - 6^ + 120^] 

T b f = \[6 H2 + ^H 2 d 

7& 6 = ±[6 + 0^ + 30/^ + 26 ] 

7*3 = 1[6 - 6>qq - 0^] 

= + 40# 2 - 20^ - 6q W + 46»^ 2 g^] 
Tti = ±[15 - 180 jj 2 g + 366 H , 2 - d^Q - 129^^ - 166 + 480^^] 
T u = ^[6 + 36 H2 + 26>£ 2 g + 0^] 
T t b = \[0hi + ^id] 
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TABLE II 



z n - ^[11 + - 80^ - 40£g + 8^g^ + 

20* - 80 Hl ^ - 20 z - 0^ z - 80^ - 40^ £ + 40^] 
Z T2 = | [-3 + 60* - 60 £ - 12fl Hl£ + 60^ - 3%^ + 120^ iZ ] 
Z rr = I [-16 + +60 #1 - #i ^ - 30* + 40^ + 40 M - 29 6lL + 80„ iZ ] 

^63 = | [6 - 29 1) + 4% - 9^ - 29 q - 40 6q - 6qq 

Zb2 = | [-3 + 69 Hl - 69 q - 129 Hi q + 69^ - 36^ + 120^£ i( =>] 

Z bl = ^[-21 + 100^ - 80^ - 4B §6 + 240^ + 180* - 240^ 

Z b6 = | [-24 + 60^ - #i6 - 30* + 40^ + 120 Q5 - 20 Q#i + 80 HiQ ] 
Z bt = ^[20^ 2 - 0* - 6fjjj 2 ] 

Z t 3 = | [6 - 20g + 40^ - 0g£ - 20^ - 40^g - O^fj] 

Z t2 = | [-3 + 60* - 60q - 120^ Q + 60^ - 30^^ + 120^ 2Q ] 

-329 - 169 M + 489 M h 2 - 486^ - 160^ Q ] 

z tb = i[20* - 9 Hl - 0_5^J 

= \[-2A9 Hi( j + 60^ 2 - 0^ - 30* + 40 H2& + 120 Q ^ - 20^ + 80 ff2Q ] 
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TABLE III 
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330.1,311.7 


330.1,311.7 


330.0,311.7 


330.0,311.6 


328.8,310.3 


£1,2 


321.0 


321.0 


320.9 


321.0 


323.5 




630.7 


612.0 


588.1 
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669.0 
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114.1,630.6 
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TABLE IV 

m t = 175 , A Q = 400 , m Q = 300 , tan/3 = 10 , fi(M z ) > 
mi/g 700 500 300 100 75 



M GUT 1.705 1.915 2.290 3.145 3.182 

(lO w GeV) 

a GUT .04030 .04066 .04120 .04215 .04222 

127.9 127.9 127.9 127.9 127.9 

sin 2 9 w .2311 .2311 .2311 .2311 .2311 

a 3 .12669 .12795 .12983 .13228 .13110 



M t 177.0 177.0 177.0 177.0 177.0 



g 1552.7 1140.1 714.6 267.8 207.1 

Xi 291.5 204.3 118.6 34.7 24.4 

X° 2 527.9 371.3 215.1 62.2 44.5 

X° 3 752.0 612.0 450.5 259.5 234.7 

Xl -733.9 -595.8 -434.8 -243.8 -219.5 

Xt 751.4 611.6 450.6 261.9 237.9 

X| 527.7 371.1 214.7 60/T 42.5 



ti,f 2 1347.6,1077.8 1022.5,780.2 698.2,470.9 388.3,167.5 353.0,144.7 

61,62 1380.8,1311.8 1039.9,975.9 695.6,634.7 370.2,305.3 337.8,270.6 

fi,f 2 548.5,390.1 446.6,346.7 362.5,313.7 311.9,294.8 308.8,293.8 
v T 542.2 438.4 351.2 296.9 293.6 



«i,2,ui )2 1444.7,1394.9 1083.2,1049.7 718.1,701.0 372.1,371.1 337.6,338.4 
di,2,d$' 2 1446.6,1389.7 1085.9,1046.7 722.1,700.3 380.0,373.4 346.3,341.1 

e li2 ,e c 12 550.0,399.3 447.3,355.5 361.8,322.8 309.1,304.7 305.9,303.7 
i/1,2 ' 544.9 440.9 353.7 299.3 295.9 



A 1341.8 1036.5 718.7 414.6 383.9 

h ,H 119.4,1341.9 118.1,1036.6 115.9,718.8 108.7,414.7 105.9,384.0 
if* 1344.0 1039.3 722.8 421.9 391.8 
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TABLE V 

m t = 175 , tan(3 = 10 , m Q = 300 , m 1/2 = 200 , fi(M z ) > 
Ao= 800 -800 300 -300 



Mqut 


2.574 


2.657 


2.645 


2.690 


2.676 


/1 nlG/^f T r\ 

(10 °GeV ) 












a GUT 


.04155 


n a "i o p* 

.04165 


.04163 


.04169 


i~\ A 1 f 'I 

.04167 


— 1 
"em 


1 r»PT r\ 

127.9 


127.9 


127.9 


127.9 


127.9 




.2311 


.2311 


.2311 


.2311 


.2311 




.13031 


.13183 


.13148 


.13235 


.13206 




1 ( {.[} 


i 77 n 


1 ( ( .U 


i 77 n 
1 f f .U 


1 77 n 


9 


493.0 


501.1 


496.7 


500.3 


498.6 


X°i 


77.3 


75.9 


76.0 


73.7 


74.6 


■CO 

x 2 


1/11 o 

141.5 


134. U 


1 OP. 1 


1 07 1 

Iz / .1 


1 on c 
130. o 


~0 

Xs 


4^7 4 


ouo.z 


OOO.O 


971 "X 


9Q1 9 

zy i.z 


Xt 


-447.5 


-282.1 


-320.5 


-240.5 


-265.7 


xl 


458.4 


305.6 


339.2 


271.4 


291.5 


xl 


141.6 


132.7 


135.3 


124.6 


128.9 


f i f n 


522 6 1 82 5 


514 Q 413 8 


541 6 332 8 


538 405 5 


543 4 378 5 


h,b 2 


520.8,433.0 


529.8,484.2 


526.8,470.9 


530.5,490.2 


529.2,484.1 


T~1,T 2 


QQi o oo/i n 
ooi.y,zy4.u 


ozo.U,oUz.o 


qqi 7 qoq rc 
ool. 1 ,oUo.O 


QQfi 1 QDfi s 

ooU.l,oUo.o 


QQfi O QDfi Q 




316 2 


317 7 


319 


319 5 


319 6 


^1,2,^,2 


536.1,527.2 


540.9,532.0 


538.2,529.2 


540.4,531.6 


539.4,530.5 


di t 2,di,2 


541.6,527.8 


546.3,532.6 


543.6,529.9 


545.8,532.2 


544.8,531.2 


ei,2,e$ >2 


330.1,311.5 


330.4,311.6 


330.0,311.7 


330.1,311.7 


330.0,311.7 


#1,2 


321.1 


321.4 


321.0 


321.1 


321.0 




661.3 


449.1 


537.2 


447.5 


480.5 


h ,H 


119.6,660.8 


105.7,449.5 


112.7,537.4 


108.2,447.8 


110.1,480.8 




665.8 


455.7 


542.7 


454.1 


486.6 
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TABLE VI 



m t 


= 110, tanp 


in A ^ 
= 1U , AO = 


400 , mi/2 = 


200 , fi(M z ) > 




8nn 

OUU 


uuu 


Ann 

1UU 


200 


100 


IVlQUT 


9 R87 
z.uo / 


9 RSI 
z.uoi 


Z.UIJU 


2.581 


2.482 


(10 w GeV) 












a GUT 


.04125 


.04138 


.04153 


.04169 


.04176 


"em 


127.9 


127.9 


127.9 


1 07 n 

iz r .y 


1 07 n 

iz ( .y 


sin t>w 


oqi 1 
.zoll 


.Zol 1 


.Zoll 


OQ 1 1 

.zoll 


OQ 1 1 
.ZOl 1 




.loU/O 


1 QAS7 


. lollU 


1 Q 1 A O 

.lol4z 


.lol4o 


M t 


177.0 


177.0 


177.0 


177.0 


177.0 


9 


510.9 


504.7 


A no A 

498.4 


/lO/l 

4y4.z 


/in/i q 
4y4.o 


X°i 


77.4 


77.1 


76.6 


7£? 

f o.z 


IP. n 


~o 

X2 


1 Q7 
lO ( .Z 


lot . 1 


1 Q7 7 


137.6 


137.5 


X°3 


340.1 


351.6 


357.3 


360.6 


361.5 


Xl 


-322.1 


-334.9 


-341.2 


Q/1 /l O 

-o44.y 


Q/1 £ O 


v c 

Al 


340 


3^9 9 


3^8 


361.4 


362.3 


X°2 


136.4 


137.0 


137.0 


1 37 n 






762.2,522.9 


652.9,432.2 


569.7,349.3 


517.5,281.9 


504.0,260.9 


bi,b 2 


879.9,734.4 


719.6,611.1 


581.4,507.1 


A 80 8 A Q1 Q 

4oz. 0,401 .y 


404.0,41U. 1 


T1 To 


805 2 793 5 


611 8 596 2 


422 6 399 7 


247.5,205.7 


178.6,113.9 




800.1 


604.9 


412.3 


99Q 9 


1 ^9 3 

lOZ.O 


~c 

U l,2) u l,2 


con 887 Q 


7Qfl 7 79c; R 


KQQ ^8^ ^ 

oyo.z,ooo.o 


494 8 484 6 


466 5 455 2 


"1,2, "1,2 


8.Q3 ^ 8.8.8 9 


70,4 r 79R q 




500.7,485.2 


472.8,455.8 


ei,2,e$ >2 


807.6,803.0 


612.8,604.9 


421.9,408.4 


244.0,217.6 


173.0,132.2 


£1,2 


804.1 


608.1 


414.9 


231.6 


155.0 




910.8 


756.5 


618.1 


514.5 


485.0 


h , H 


112.3,911.0 


112.4,756.7 


113.2,618.2 


114.4,514.6 


114.8,485.0 




914.1 


760.4 


622.9 


520.3 


491.1 
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TABLE VII 



tan(3 = 20, Ao 


= zUU , m = 


400 , mi/2 = 


300 , fi(M z ) > 


m t 


175 


1 8n 

lOU 


1 8^ 


190 


195 


Mqut 


2.337 


9 ZL1 ZL 


9 ZLQ^ 


2.581 


2.670 


(10 w GeV) 












a GUT 


.04117 


.04118 


.04119 


.04120 


.04120 


^ -1 

a em 


iz ( .y 


127.9 


127.9 


iz i .y 


i 07 n 

iz i .y 


sin 2 6w 


oq 1 1 n 


.Zoliy4 


.ZoU ( O 


OQf|/?1 

.ZoUol 


.ZoU44 




.loUUo 


.loUoo 


1 Q 1 1 /I 


.lol 1 1 


.lozzO 


M t 


177.0 


181.7 


186.4 


191.1 


195.8 


9 


717 q 


719.0 


720.7 


799 £ 
1 ZZ.O 


79/1 A 
1 Z4.4 


Xi 


iion 

iio.y 


119.0 


119.1 


iiy.z 


11(1 Q 


X°2 


213.2 


91 ^ n 


91 /I 
Z10.4 


217.5 


218.5 


X°3 


402.0 


422.0 


441.3 


459.9 


477.9 


~o 

Xi 


-oo4.U 


-406.2 


-427.2 


A A 7 9 
-44 ( .Z 


-400. 


xl 


402.6 


499 7 


449 1 


460.8 


478.8 


X2 


91 9 8 


214.7 


216.1 


917^ 


91 8 3 




720.3,526.9 


717.7,522.9 


715.2,519.5 


712.9,516.8 


710.7,515.1 


h h 


797 D KKK 1 
I Z I .U,OOD. 1 


728.9,663.1 


730.8,661.3 


739 7ffiQ8 

i oz. i ,ooy .o 


73/1 fi^S 7 
1 O4.0,00o. ( 


fi,r 2 


444.4,392.8 


444 9 392 2 


445 5 391 6 


446.0,391.1 


446.5,390.5 




4^9 9 


432.3 


432.4 


439 ^ 


439 ^ 


"1,2; "1,2 


761 3 745 7 


7R9 7 747 9 

( OZ. 1,141./ 


7R4 3 748 7 


765 9 750 3 


767 5 752 


di,2,di,2 


765.1,745.2 


7fifi fi 74fi fi 


7fi8 1 748 9 


769.7,749.8 


771.3,751.5 


ei,2,e$ >2 


447.3,417.0 


447.3,417.0 


447.4,417.0 


447.4,417.1 


447.5,417.1 


#1,2 


440.7 


440.8 


440.8 


440.9 


441.0 




795.6 


820.6 


843.0 


862.8 


879.8 


h ,H 


115.0,795.6 


116.9,820.6 


118.8,843.0 


120.7,862.7 


122.6,879.7 




799.3 


824.1 


846.4 


866.1 


883.1 
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TABLE VIII 


m t = 175, tanf3 = 10, A Q 


= ZOU, 77l = 


9nn «n i en i,(i\/f \ \ n 
ZUU, — loU, /i^M^J > U 




Case [a] 


Pqco [hi 


oase [cj 




1-loop 


2-loop 


Complete 2-loop 


M G ut 


2.1881 


9 887K 
Z.oo 1 O 


9 87RR 


(10 w GeV) 








OLGUT 


.04127 


.04201 


.04202 


„ -1 

"em 


iz ( .y 


127.9 


127.9 


• 2/1 
Sin Uyy 


.zolUo 


9Q 1 1 n 


9Q1 1 n 


«3 


.11 ( tw 


1 Q9Q A 




M t 


181.0 


177.0 


177.0 


9 


one /i 


381.4 


382.6 


~o 

Xi 


Dy.Z 


55.0 


54.4 


X°2 


109.0 


yo.o 


Ofi 7 

yo. i 


xl 


302.8 


304.1 


279.6 


Xi 


OS/1 1 

-zo4. 1 


-287.7 


-260.1 


Xi 


304.0 


30^ 


980 Q 


~c 
X2 


1UO.U 


97.4 


95.3 


tj 2 


443.6,247.0 


442.0,235.2 


440.2,234.7 


h h 


/im 7 qc;7 9 

4U1. 1 ,oO I ./ 


395.4,352.8 


394.9,352.6 


Tl,f 2 


235.3,203.6 


232 2 201 6 


231 3 202 7 




91 fi fi 
Z10.0 


212.6 


212.5 


u 'l,2)"'l i 2 


410 5 401 4 


/inn i ^Q/i i 


/inn i ^Q/i 1 


^1,2,^1,2 


417.8,402.4 


407 c 3Q C 7 


4fl7 c 30,^ 7 


ei,2,e^ 2 


231.6,212.6 


227.7,211.6 


227.5,211.8 


#1,2 


218.1 


214.2 


214.1 




412.1 


421.0 


394.5 


h ,H 


113.0,412.2 


110.7,421.1 


110.5,394.7 


H± 


419.5 


428.1 


402.1 
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